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Abstract 

The Navier-Stokes equations, are understood as the result of the low-order 
expansion in powers of dimensionless rate of strain ijij = ToSij, where tq 
is the microscopic relaxation time of a close-to- thermodynamic equilibrium 
fluid. In strongly sheared non-equilibrium fluids where \riij\ > 1, the hydrody- 
namic description breaks down. According to Bogolubov's conjecture, strongly 
non-equlibrium systems are characterized by an hierarchy of relaxation times 
corresponding to various stages of the relaxation process. A " hydro- kinetic" 
equation with the relaxation time involving both molecular and hydrodynamic 
components proposed in this paper, reflects qualitative aspects of Bogolubov's 
hierarchy. It is shown that, applied to wall flows, this equation leads to qual- 
itatively correct results in an extremely wide range of parameter ry-variation. 
Among other features, it predicts the onset of slip velocity at the wall as an 
instability of the corresponding hydrodynamic approximation. 
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1 Introduction 



Strongly sheared fluids, in which the usual Newtonian hydrodynamic description 
breaks down, are commonly encoimtcrcd in biology, chemical engineering, micro- 
machinery design [l]-[2]. Extensive efforts that largely relied upon physical intuition 
and qualitative considerations to incorporate corrections at the hydrodynamic level 
of description have been made during the years and have achieved various successes. 
However, these attempts generally fail when the non-linearity of a fluid is strong. 
Furthermore, theoretical understanding of these highly non-linear systems remains a 
major challenge. 

The Navier-Stokes equations, which can be derived from the Boltzmann kinetic 
equation as a result of a low -order truncation of an infinite expansion in powers of 
dimensionless length-scale Kn , have been extremely successful in describing Newto- 
nian fluid flows [3]- [7]. The parameter Kn is defined as a ratio between the so called 
relaxation time tq associated with molecular collisions and the time of the molecular 
convection, namely 

Kn = Tq/{L/Cs) = TqCs/L 

where L is a characteristic length scale a flow inhomogeneity, and is the sound 
speed (i.e., an average speed of the molecules). According to kinetic theory, parameter 
A = tqCs represents the so called mean free path in a system and tq is the time-scale for 
the system to relax to its local thermodynamic equilibrium. Since the ratio L/cs is a 
characteristic time-scale of deviation from equilibrium due to density (concentration) 
perturbations, the Knudsen number Kn is a measure of departure of a fluid system 
from thermodynamic equilibrium in [4] -[6]. 

For the purpose of understanding the effects of non-linearity, it is desirable to 
re-express Kn in terms of the rate of strain (velocity gradient) in the fluid. Using the 
estimation for the hydrodynamic (macroscopic) length scale, 

L^U/\S\ 

we can also write Kn in an alternative form, 

Kn = T^\S\/M, = r)o/M, (1) 

In the above, U represents the characteristic velocity, and \S\ is the magnitude of the 
rate of strain tensor Sij of the flow fleld (l^*! = ^2SijSij). Ma = U/cg is the Mach 
number. The rate of strain tensor is commonly deflned as, 

_ 1 duj duj 
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The dimensionless rate of strain, t]q = tq\S\ (in what follows we, instead of ?7o, will 
use the parameter t] = tS where r is a physically relevant relaxation time) can 
be viewed as a measure of degree of inhomogeneity and shear in the flow. It is 
equivalent to a parameter rj ^ W = Utq/L commonly used in hydrodynamics of 
polymer solutions. Eqn.(0) indicates that Kn can be large if t] is large, specially 
for low Mach number flows. In many situations rj is not a small parameter and the 
Newtonian fluid based description breaks down. For example, capillary flows or blood 
flows through small vessels, rarified gases and granular flows cannot be quantitatively 
described as Newtonian fluids. A finite rj can either be a result of a large relaxation 
time associated with the intrinsic fluid property (such as in some polymer solutions), 
or a result of strong spatial variations like turbulent and micro or nano-flows. This 
latter effect is particularly important at a solid wall where the velocity gradient is 
often quite significant. It is known that the no-slip boundary condition is only valid 
in the limit of vanishing Kn. Experimental data indicate that the no-slip condition 
breaks down when Kn > 10~^ [8] and the Navier-Stokes description itself becomes 
invalid at Kn > 0.1. For example, the experimentally observed velocity profile in a 
simple granular Couette flow does not resemble the familiar linear variation of velocity 
predicted from the Navier-Stokes equation with no-slip boundary conditions [9]. The 
wall slip is an indication of an extremely strong local rate of strain. 

The Navier-Stokes equations can be perceived as a momentum conservation law 
with a linear (Newtonian fluid) stress-strain relation [3], [5]- [7], i.e. the deviatoric 
part of the stress tensor takes on the following form, 

= (v'iv'j) - ^{V^)Sij ^ uSij, i,j = x,y,z (2) 

where coefficient u is the kinematic viscosity. In the above, v' (= v — u) is the relative 
velocity between the velocity of molecules (v) and their locally averaged velocity (u) 
namely the fluid velocity. It is expected that v' and u represent, respectively, the fast 
(kinetic) and the slow (hydro dynamic) velocity fields. In an unstrained flow where 
the rate of strain is equal to zero {Sij = 0) , the relation (0) is simply interpreted as 
the first term of an expansion in powers of small rate of strain. Therefore, to describe 
rheological or micro-flows with high rate of strains, such a linear approximation must 
be modified to include non-linear effects. However, this task is highly non-trivial if 
not impossible. 

The hydrodynamic approximations can be derived from kinetic equations for the 
distribution function with the intermolecular interactions accounted for through the 
so called collision integral [3], [5]- [7]. One can formally write the Bogolubov chain 
of equations for the distribution function and, in addition to the high powers of the 
dimensionless rate-of-strain rj, generate an infinite number of equations for the multi- 
particle contributions to collision integral [3], [10]. For a strongly sheared flows where 
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the expansion parameter t] is of order unity, this chain cannot be truncated. An 
additional difficulty is that the consistent expansion includes the high-order spatial 
derivatives , which unfortunately means that, even if we were able to develop the 
procedure, the resulting un-truncated infinite- order hydrodynamic equation requiring 
an infinite set of initial and boundary conditions would be useless. 

In this paper, based on Bogolubov's concept of the hierarchy of relaxation times, 
we propose a compact representation of the infinite- order hydrodynamic formulation 
in terms of a simple close-form hybrid ( "hydro-kinetic" ) equation. The power of the 
approach is demonstrated on a classical case of stronly sheared non- Newtonian fluids. 
It is shown that in our "hydro-kinetic" approach, the formation of a slip velocity 
at a solid wall and the simultaneous flattening of velocity distribution in a bulk is 
a consequence of an "instability" of the corresponding hydrodynamic equation with 
the no-slip boundary conditions. This "instability" is the result of the non-universal 
details of the flow such as the local values of dimensionless rate of strain r]. 

2 Basic formulation 

Following the standard Boltzmann kinetic theory, we introduce a density distribution 
function /(x, v, t) in the so called single particle phase space (x, v) (= F). A formally 
exact kinetic equation, involving an unspecified collision integral C, can be given [11] 

dJ + v,^f = C (3) 

Generally speaking, the detailed expression for C involves multi-particle distribution 
(correlation) functions /^"^(ri, . . . , F^, t) (with n > 1). This results in the famous 
Bogolubov chain of infinite number of coupled equations [3], [10] . This chain cannot 
be closed when fluid density is not small. On the other hand, the collision integral 
C can be modeled in a relatively simple form in the case of a rarified gas where only 
the binary collisions are important [ 11] 

C = l w'if'f[ - f h) dT,drdT\ = I Urei{ff[ - f fi)dad'p (4) 

In this model, C depends only on the single particle distributions. Here C dT is the 
rate of change of the number of molecules in the phase volume dT = d^xd^v, Urei is 
the relative velocity of colliding molecules and da and p are the differential scattering 
cross-section and momentum of the molecules, respectively. The state variable T 
describes all degrees of freedom of a molecule and w' stands for the probability of a 
transition of two molecules initially in states T and Fi to states F' and Fi' as a result 
of collision. The kinetic equation Q together with the specific collision integral (j3]) 
forms the celebrated Boltzmann equation. 
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It is well known that Boltzmann equation admits an H-theorem in that the system 
monotonically decays (relaxes) to its thermodynamic equilibrium. The corresponding 
local thermodynamic equilibrium distribution function, f'^'^, is determined from the 
solution C{f^'^) = 0. If deviation from the equilibrium is small, we can write / = 

+ 6f and 

C = -^+ ureif'f[dad'p ^ — (5) 

To J Tq 

with ^ 

To J 

Equation © is often referred to as the "BGK" (mean-field) approximation [12] which 
is a natural reflection of the Boltzmann H-theorem. Furthermore, when perturbation 
from equilibrium is weak, it indicates that the relaxation to equilibrium is realized for 
each of the distribution functions individually having a common relaxation time. At 
this point, it is important to make the following clarification: Even though as shown 
above that was deduced from the Boltzmann binary collision integral model, its 
applicability can be argued to go beyond the low density limit. Indeed, the model 
process is consistent to the more general principle of the Second law of thermody- 
namics: A perturbed fluid system always monotonically relax to thermal equilibrium, 
regardless whether the system has low or high density. Furthermore, in accord with 
Bogolubov [10] the general process of return-to-local equilibrium is true even when 
the deviation from equilibrium is not small. 

In the classical kinetic theory, the relaxation time tq (~ const) represents a char- 
acteristic time of the relaxation process in a weakly non-homogeneous fluid and the 
smaller value of tq, the faster the process of return to equilibrium. The formal expan- 
sion of kinetic equations in powers of dimensionless rate of strain ( Chapman- Enskog 
(CE) expansion [3], [6]) developed many years ago, leads to hydrodynamic equations 
for the macroscopic ("slow") variables. When applied to ©, the first order trancation 
of the formally infinite expansion gives the Navier-Stokes equations with kinematic 
viscosity: 

z/Q = tqO. (6) 

Development of the expansion based on the full Boltzmann equation is an extremely 
difficult task. On the other hand, the simplified Boltzmann-BGK equation (5) for the 
single-particle distribution function allows extension of the Chapman-Enskog expan- 
sion to include higher powers of dimensionless parameter rj. Considering the higher- 
order terms of the Chapman-Enskog expansion, a simple scalar relaxation term in (5) 
combined with advection contribution, is expected to generate an infinite number of 
anisotropic contributions as contracted products of Sij. Indeed, by expanding equa- 
tion (5) up to the second order in the Chapman-Enskog series, we can explicitly show 
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that the deviatoric part of the momentum stress tensor takes on the following form 
[13] 



2T9S,j-2Te{dt + u-V){TSij) 

—4:T^9[SikSkj — -jSijSkiSki] + 2T'^9[SikQ,kj + Sjk^ki] 



(7) 



where the vorticity tensor is defined as, 




1 p dui duj 



2 dxj dxi 



Note the first 0{r]) term in (7), resulting from the first order Chapman- Enskog (CE) 
expansion, corresponds to the Navier-Stokes equations for Newtonian fluid, while the 
non-linear corrections to the Navier-Stokes hydrodynamics are generated in the next 
(0(77^)) order of the CE expansion. It is important to further point out that the 
memory effects, which appeared in the hydrodynamic approximation (7) as a result 
of the second-order CE expansion, are contained in a simple equation (5) which can 
be regarded as a generating equation for hydrodynamic models of an arbitrary- order 
non-linearity . 

The hydrodynamic approximation (6), (7) has been derived from the equation 
(5) for the single-particle distribution function valid for a weakly non-equilibrium 
fluid (small rf). In strongly sheared fluids both the assumption of local equilibrium 
and the low-order trancation of the Bogolubov chain (single-particle coUisions) fail 
and the accurate resummation of the series is impossible. Indeed, even in the low- 
density fluids, the strong shear {rj >1) introduces local fluxes facilitating long-range 
correlations between particles, invalidating the fluid description in terms of the single- 
particle distrubution functions. To account for these effects, we have to modify kinetic 
equation (5) in accord with some general ideas about the relaxation processes. 
Our goal is to reformulate the kinetic equation (5), by modifying the relaxation time 
r and come up with the effective kinetic equation qualitatively accounting for the 
effects of the neglected high-order contributions to the Bogolubov chain dominating 
the dynamics far from equilibrium. 

In his seminal 1946 book [10] Bogolubov proposed the hierarchy of the time-scales 
that describe relaxation to equilibrium for a system initially far from equilibrium. 
According to his picture, these initially strong deviations from equilibrium rapidly 
decrease, thus allowing an accurate representation of the entire collision integral in 
terms of the single-particle distribution functions. This dramatically simplified repre- 
sentation is sufficient for an accurate description of the later, much slower, process 
of relaxation to thermodynamic equilibrium. 
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To make this plausable assumption operational, we have to represent a Bogolubov 
hierarchy of relaxation times in terms of observable dynamical variables characterizing 
the degree of deviation from equilibrium. Since the most natural parameter governing 
the dynamics far from equilibrium is the rate of strain \Sij\, we propose that both 
the close-to-equilibrium relaxation time tq and inverse rate-of-strain l/li^jjl define the 
Bogolubov hierarchy of the relaxation times. The simplest Galileo invariant relaxation 
model that is compatible with the above physical considerations is: 



where we define in accord with (2), S^J^ ~ ^ and rjne = to\S"''^\ {i ^ j). In what 
follows, wherever it cannot lead to misunderstandings, we will often omitt the suffix 
ne. The new collision integral (i.e., (jSj) with tq replaced by r in (jS))) now describes a 
relaxation process that involves a rate of strain-dependent relaxation time. 

The proposed hydro-kinetic model (^-^) is chosen to reflect some of the princi- 
ple elements in the Bogolubov hierarchy. That is, far from equilibrium where the rate 
of strain is large, the essential time-scale t is dominated by \S\~^ which corresponds to 
a rapid first stage of the relaxation process. Later on, when the rate of strain becomes 
small, the relaxation process is governed by a close-to-equilibrium relaxation time tq, 
as used in the conventional BGK equation (5) leading to the Navier-Stokes formula- 
tion. It is clear from relation (7) that even though, by restricting our model to the 
scalar relaxation times in (jSJ, the anisotropic contributions to the stress do appear in 
the hydrodynamic description which is a result of the Chapman- Enskog expansion. 
Since the rate of strain is a property of the flow, the model (8) which includes both 
molecular and hydrodynamic features can be called a hybrid hydro-kinetic approach 
to non-linear fluids. It is shown below that even with such a specific form in (jHJ, the 
hydro-kinetic equation is capable of producing some quite non-trivial but physically 
sensible results at the hydrodynamic level. 

To sumamrise the main points of this section: it is analytically impossible and 
practically useless to attempt to describe the strongly non-linear (non-equilibrium) 
fiow physics at the hydrodynamic level. As indicated above, not only the resulting 
"differential equation" does not have a finite form, it also requires an infinite number 
of boundary conditions. The reason for this fundamental difficulty is that the expan- 
sion in powers of parameter t] becomes invalid when r] is not small. Thus, to deal 
with strongly sheared and time-dependent non-linear fiows, it is desirable to use the 
simple "hydro- kinetic" description Q-dHl)- Clearly, this hybrid representation has a 
finite form, while it formally contains all the terms in the infinite expansion for the 
hydrodynamic level. As argued above, the hydro-kinetic formulation is applicable to 
both large and small rj, corresponding to large and small deviations from equilibrium. 
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3 Wall flows of non-Newtonian fluids 



To illustrate the benefits of the hybrid representation of (jSI)-®, let us first consider a 
laminar unidirectional flow in a channel between two plates separated by distance 2H 
and driven by a constant gravity force g. In a steady state, the Navier-Stokes equation 
for the channel flow can be derived from (5), (8) in the lowest order of the expansion 
in powers of dimensionless rate of strain rj = T\Sij\. Repeating the procedure leading 
to (7) and restricting the expansion by the first term gives the Navier-Stokes equation 
having a "renormalized" viscosity corresponding to (jHl) , 

dtu — g = dyvdyU (9) 

where 

v = re = p^/^l + -iTi\S\^ (10) 

for the special unidirectional situation, 15*1 = \du/dy\. In the above, we have chosen 
u to be the streamwise velocity component, while coordinates x and y are along the 
streamwise and normal directions of the channel, respectively. 

One important point must be mentioned: the expansion in powers of r/o = tqS 
corresponds to the classic CE expansion. The equation (9) has been derived by 
expanding in powers of = t\S\ = 7]^/ \Jl + 7779, which means that even the low- 
order perturbation theory in powers of the " dressed " parameter r] corresponds to an 
infinite series in powers of "bare" parameter rj^. Thus, it is extremely interesting to 
assess the accuracy of the derived hydrodynamic approximation (9)-(10). 

Subject to no-slip boundary conditions, the exact steady-state analytic solution 
for Q is given by: 

where (3 = '-/{rogH/uo)^ is a dimensionless parameter which can be either positive 
or negative, for 7 can in principle have either positive or negative signs. One can 
immediately see that this particular flow solution breaks down for (3 > 1. A direct 
indication of this is that the steady state Navier-Stokes equation © is no longer 
valid to describe such a flow, and we have to return back to the full hydro-kinetic 
representation (jSl)-®- It is interesting that the equation (9)-(10) allows an unsteady 
singular solution: 

u = gt + Uo(f){y) (12) 

where 0(?/) = 1 for \y\ < H and (f){y) = on the walls where \y\ = H. The transition 
between the two (no-slip (11) and slip (12)) solutions will be demonstrated below. 
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In the rest of the paper, we present solution of the full hydro-kinetic system for the 
channel flow for the entire range of parameter variation of 7 with initial velocity profile 
u{y,t = 0) = 0. For this purpose, equations (jSI)-® have been numerically solved 
using a Lattice Boltzmann (LB) algorithm [14]. On each time step the relaxation 
time in (8) was calculated with the non-equilibrium rate of strain S^J^, defined as: 

SJ^; = - f dv,dv,{v, - u,){v, - u,)U - D (13) 

with i ^ j and u standing for the local value of the mean velocity. It is clear from 
this definition that in thermodynamic equilibrium, the rate-of-strain S'fJ = and, 
according to (7), not far fom equilibrium Sij ~ S'^-^. The computationally effective 
and widely tested "bounce-back" collision process giving, on the hydrodynamic level 
of description, rise to the no-slip boundary condition in the Kn limit, was imposed 
on a solid wall. According to the "bounce back" algorithm, the momentum of the 
"molecule" colliding with a solid surface changes according to the rule: p ^ — P- 

When parameter /? = , the familiar steady state parabolic solution u{%j) = 
^^(1 — (-^^) was readily derived. Figs.l present the analytical (i.e., (fTTjl ) and the 
numerical (i.e., 0-®) solutions of the velocity profiles in the plane channel flow 
for, respectively (3 = 0.4. As we can see, for this value, the difference between the 
solutions of (fTTj) and simulations of the full hydro- kinetic model (©-(0)) is negligible. 
This means that in this regime, the lowest order trancation of the Chapman-Enskog 
expansion in powers of the "dressed" parameter t] (7) is extremely accurate. The 
same conclusion was shown to hold for all (3 < 0.5. 

The numerical results for (3 > 0.51 revealed an interesting instability theoretically 
predicted , for (3 > 1. The results of very accurate numerical simulation (960 cells 
across the channel) are presented on Figs. 2-3. 

We can notice qualitatively new features not captured by the steady-state hy- 
drodynamic approximation: initially, formation of aly narrow wall-boundary layer, 
accompanied by a strong flattening of the velocity profile in the bulk can be observed. 

Later in evolution, the boundary layer becomes unstable with formation of the slip 
velocity at the wall. The flow accelerates, eventually becoming a free -falling plug 
flow, predicted by equation (12). All these phenomena have been experimentally 
observed in the flows of rarified gases and granular materials. 

A clarification of the set up is in order: The simulations were performed on an ef- 
fectively infinitely long (periodic boundary conditions along the streamwise direction) 
channel, and the flow was driven by an externally imposed constant gravity. This set 
up differs quite substantially from a pressure-gradient-driven flow of a nonlinear fluid 
where a steady state can be achieved by formation of a non-constant (x-dependent) 
streamwise pressure gradient. Unlike pressure, gravity is not a dynamical variable 
and hence the flow lacks the mechanism for establishing a force balance needed to 
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Figure 1: Steady state velocity profile: comparison of hydrodynamic solution Eq. 11 
(x) with the LBM simulation o. (3 = OA. 

achieve a steady state. This can be associated with the experimentally observed in- 
ability of the gravity-driven granular flows in ducts to reach steady velocity proflles 
[14]. In accord with this theory, the steady velocity proflle similar to those shown 
on Figs. 2-3 can also be observed in a gravity- driven finite-length-pipe or channel 
flows . In this case we expect the velocity distribution to vary with the length of the 
pipe/channel. 
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Figure 2: Short-time evolution of unstable velocity profile (3 = 0.51. The values of 
dimensionless time in arbitrary units: from bottom to top: T=0; 2; 4; 6. 

4 Conclusion 

It has recently been shown that the Lattice Boltmann ( BGK ) equation with the 
effective strain-dependent relaxation time can be used for accurate description of 
high Reynolds number turbulent flows in complex geometries [16]. In this work, this 
concept has been generalized to flows of strongly non-linear fluids. Although the 
simple relaxation model (jH} was proposed here on a qualitative basis, it has shown to 
be capable of producing non-trivial predictions for flows involving strong non-linearity. 
To the best of our knowledge, this hybrid ("hydro- kinetic") model is the flrst attempt 
of incorporating the principle elements of the Bogolubov conjecture about inflnite 
hierarchy of relaxation times. The most interesting result of application of this model 
is the appearance of the slip velocity on the wall as a result of dynamic transition 
driven by increasing rate of strain. Since this transition depends on the wall geometry, 
it cannot be universal. Thus, to predict this extremely important effect, the model (3)- 
(8) does not require empirical, externally imposed boundary conditions. The classical 
incompressible hydrodynamics relies upon one externally determined parameter, the 
viscosity coefficient which can be obtained either theoretically (sometimes) or from 
experimental data. The hydro-kinetic approach proposed in this paper needs a single 
additional parameter 7 describing physical properties of a strongly non -linear fluid 
, which can readily be established from a low Reynolds number flow in a capillar by 
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Figure 3: Long time velocity profile evolution (3 = 0.51. The values of dimensionless 
time in arbitrary units: from bottom to top: T=0; 8; 16;32; 80. 

comparing the measured velocity profile with the theoretical prediction (11). Further 
applications of the model (3)-(8) to the separated highly non-linear flows, will show 
how far one can reach using this simple approach. 
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